Electronic structure of fullerene-like cages and finite nanotubes of aluminum nitride 
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We report density functional study of alternate fullerene-like cage structures and finite closed, 
capped single-wall nanotubes of aluminum nitride. The cages and nanotubes studied are modeled 
as AI24N24, AI28N28, AI32N32, AI36N36, AI48N48, and AI96N96. The structure optimization and 
calculation of the electronic structure, vertical ionization potential, and the electron affinity are 
performed at the all electron level by the analytic Slater-Roothaan method, using polarized Gaussian 
basis set of double zeta quality. All structures are energetically stable with binding energy of about 
10-11 eV per A1N pair. For the larger AI96N96, the fullerene like cage is energetically less favorable 
than the two-shell cluster that has AI24N24 as an inner shell. The vertical ionization potential and 
the electronic affinity are in the range 6.7-6.9 eV and 1.5-2.0 eV, respectively. The binding energy 
show systematic increase with increase in the length of (4,4) nanotube. The energy band gap, 
determined using the ASCF method show that these structures are characterized by a fairly large 
band gap about 4 — 5 eV, which is however smaller than the gap for the corresponding boron nitride 
structures. 
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I. INTRODUCTION 

Aluminum nitride (A1N) belongs to the group III-V 
family of semiconductors and possesses number of prop- 
erties such as low thermal expansion (close to that of 
silicon), high thermal conductivity, resistance to chemi- 
cals and gases normally used in semiconductor process- 
ing, and good dielectric properties, which make it po- 
tential candidate for technological applications in micro- 
electronics. Unlike another group III-V semiconductor, 
boron nitride (BN), A1N does not exist in the layered 
graphitic form, which can be rolled up to form nan- 
otubes. Nevertheless, a few reports of A1N nanotube and 
nanowire synthesis have appeared in literature in recent 
y earg> i i 2 i 3 J 4 i 5 j n ^ ne ear ii es t report^ the A1N nanotubes 

(A1N-NT) and nanoparticles were synthesized by dc-arc 
plasma method and were characterized by the transmis- 
sion electron microscope (TEM). The nanotubes were 
measured to be 500-700 nm in length with 30-200 nm 
diameter, and the nanoparticles had diameters in the 
5-200 nm size range. The field emission pattern from 



tungsten tips coated with A1N tubes was measured and 
was attributed to the tubes having open ends. Using 
the same dc-arc plasma method, helical and twisted A1N 
nanotubes were reported by the same group^ The nan- 
otubes were dispersed on graphite substrate and analyzed 
by scanning tunneling microscopy. The tubes were re- 
ported to have an average diameter of 2.2 nm and lat- 
eral dimension of about 10 nm. The interatomic dis- 
tance between two Al-Al or N-N atoms was measured 
to be 3.2 A and from I-V curves the A1N-NT were ar- 
gued to be metallic in character^ In another study, 3 
same group reported synthesis of A1N nanoparticles in 
the size range 15-18 nm and of nanowires of 500-700 nm 
in length with diameters in the 30-100 nm size range. 
In a subsequent report, these nanowires were interpreted 
to be nanotubes^ Synthesis of faceted hexagonal A1N 
nanotubes has also been reported^ These A1N-NT were 
characterized by TEM and are found to be of a few mi- 
crometers in length and 30-80 nm in diameter. Most 
tubes were found with both ends open. There are also 
a few reports of theoretical calculations £L£ Zhang and 
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Zhang 6 considered two model AI27N27 structures, one 
for an A1N nanowire and one for an A1N nanotube at 
the Hartree-Fock level using the 3-21G basis set. They 
noted that A1N bond length in finite AIN-NTs, modeled 
by the AI27N27 tubule, alternates between 1.76 A and 
1.77 A and posses large energy gap of 10.1 eV between 
the highest occupied molecular orbital (HOMO) level and 
the lowest unoccupied molecular orbital (LUMO) level. 
This rather large gap of 10.1 eV is an artifact of the 
Hartree-Fock approximation which models exchange ef- 
fects exactly but ignores correlations effects. Our cal- 
culations, on the other hand, calculate the band gap by 
the so called ASCF where the first ionization poten- 
tial is subtracted from the first electron affinity. In den- 
sity functional theory for finite systems, this is known 
to provide better approximation to the band gap than 
the significantly underestimated HOMO-LUMO eigen- 
value difference. Another theoretical calculation by Zhao 
and coworkers^ have used density functional theory in 
the local density approximation (LDA) and generalized 
gradient approximation (GGA) to study strain energy 
(energy required to curl up nanotube from graphite like 
planar sheet) and stability of selected single-walled A1N- 
NT. These calculations employed a localized numerical 
orbital basis and psuedopotential for description of ion 
cores. These authors noted HOMO-LUMO gaps of 3.67 
and 3.63 eV for (5,5) and (9,0) tubes, respectively. Sim- 
ulations at elevated temperatures at the level of LDA 
indicated stability of single walled A1N-NT at room tem- 
perature. In the present paper, we report the study of 
fullerene analogue of A1N cages and finite (4, 4) single- 
walled A1N nanotubes using analytic density functional 
theory. Unlike A1N-NT, only one study has so far ad- 
dressed A1N nitride cage structures. The study was per- 
formed using density functional theory and indicated pos- 
sible existence of AI12N12, AI24N24, and AL30N60 on the 
basis of energetics and vibrational stability, 9 The present 
study is motivated partly by these works and partly by 
reports of synthesis of boron nitride (BN) cages ^ILl^ll 
Here, we study several fullerene analogues of A1N and se- 
lected single- wall A1N nanotubes containing up to about 
200 atoms. Some of the cage structures studied in this 



work have been proposed as candidate structures for 
boron nitride (BN) cages synthesized and detected in 
mass spectrumiifliiiiiSiii Our study, performed at the all 
electron density functional level, does not preclude the 
existence of cage structures for the A1N, in agreement 
with conclusion drawn by Chang and coworkers^ The 
binding energy, the first ionization potential, and elec- 
tron affinity are calculated. The band gap calculated 
by difference of self-consistent solution (ASCF) method 
show that the cages and tubes of A1N, like their bulk 
phase, are characterized by large gap. 

II. METHODOLOGY AND COMPUTATIONAL 
DETAILS 

Our calculations are performed using the Slater- 
Roothaan (SR) method. 14 It uses Gaussian bases to 
fit the orbitals and effective one body Kohn-Sham po- 
tential of density-functional theory^ The SR method 
through robust and variational fitting is analytic and 
variational in all (orbital and fitting) basis setsi 1 ^ The 
most general functionals that it can treat so far are cer- 
tain variants 16 of the Xa functional iLHL In particular, 
it can handle different a's on different elements analyt- 
ically and variationally so that the atomized energies 
of any cluster can be recovered exactly, and all ener- 
gies are accurate through first order in any changes to 
any linear-combination-of-atomic-orbitals LCAO or fit. 
The s-type fitting bases are those scaled from the s-part 
of the orbital basis^ A package of basis sets has been 
optimized^ for use with DGauss^i We use the valence 
double-^ orbital basis set DZVP2 and the pd part of 
the (4,3;4,3) (A2) charge density fitting basis. We use 
the a values of 0.748222 and 0.767473 for Al and N, 
respectively.- 2 Using these a values with the above ba- 
sis sets, exact total energies of Al and N atoms can be 
obtained. Hence, in the present computational model, 
the atoms in the cages or tubes will have exact atomic 
energies in the dissociation limits The accuracy of this 
method using the exact atomic values of a as judged from 
the binding energies of Pople's G2 set of molecules lies 
between that of the local density approximation and that 




FIG. 1: (Color online) Two different views of optimized oc- 
tahedral O, Si, and Ss AI24N24 cages. 

of the generalized gradient approximation^ Its main ad- 
vantage is that it requires no numerical integration and 
hence within the accuracy of the model gives results 
that are accurate to machine precision. The optimiza- 
tion is performed using the Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) algorithm 2i ?i 2 ^ 2 7i 2 ?i 2 ? with forces on the 
atoms computed non-recursively using the 4-j generalized 
Gaunt coefficients. 24 

III. RESULTS AND DISCUSSION 

Before we present our results on larger A1N cages, 
we test the computational model for the simple di- 
atomic A1N molecule, for which experimental values 
of bond length and dissociation energy are available. 



FIG. 2: (Color online) Favorable T A1 28 N28 and T d A1 3 6N 36 
cages containing six square. 

TABLE I: The bond length, dissociation energy of 3 II A1N 
molecule obtained by present analytic Slater- Roothan method 
in comparison with those obtained from some selected com- 
putational models and experiment. The present value of dis- 
sociation energy does not contain zero point energy. 



Model 




Re (A) 


D (eV) 


Expt. 


ReOO 


1.79 


2.86 ±0.39 


CCSD(T)/WMR 


Re02 


1.79 


2.45 


MRCI 


ReOl 


1.83 


2.42 


MRCI 


Re03 


1.81 


2.20 


BP86/6-31G* 


Ref.34 


1.81 


2.87 


BPW91/6-31G** 


Re05 


1.80 


2.76 


Present 




1.81 


2.71 



A number of theoretical calculations at different level 
of sophistications have also been performed on this 
moleculei ? 1 ^ 2 '?^? 4 '^ The comparison of these with the 
numbers obtained in present model is given in Table [IJ 
where the present values are found to be in good agree- 
ment with the experimental values. All theoretical val- 
ues of dissociation energy are within the rather large ex- 
perimental error bar. Predicted dissociation energy in 
present model lies between the ab initio coupled clus- 
ter (CCSD) and multi reference Configuration interac- 
tion MRCI numbers on the one hand and the Becke- 
Perdew86(BP86) or Becke-Perdew-Wang91 (BPW91) 
density functional models on the other. Use of larger 6- 
311G** orbital basis and the Turbomole fitting basis used 
in our earlier workSS^ does not change the results. Bond 
length is practically identical and the binding energy de- 
creases by 0.2 eV. The basis set effect is largely canceled 
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by adjusting a to get exact atomic energies. Thus, we 
expect present model to provide a reliable description of 
A1N cages and nanotubes studied in this work. 

The optimized AI24N24 clusters of octahedral, S4 and 
58 symmetry are shown in Fig. ^ The octahedral cage 
consists of 6 octagons, 8 squares, and 8 hexagons. It is a 
round cage with two symmetry inequivalent atoms. This 
round cage has been studied as a candidate for cages 
made up of carbon^ silicon hydrides^ boron oxides?^ 
aluminum hydrides^ etc. It has also been proposed as 
structure of the B24N24 peak, found in abundance in 
a recent experimental time-of- flight, mass spectrometric 
study iii The optimized coordinates of the AI24N24 cage 
for two inequivalent atoms in A are (0.868, 2.101, 3.326) 
for Al and (3.48, 2.13, 0.916) for N, respectively. Posi- 
tions of other atoms can be determined from the symme- 
try group operations of an octahedral point group with 
symmetry axes along the coordinate axes. The cage has a 
radius of 4.11 A, about 27% larger than that of the octa- 
hedral B24N24 cage< 3 2, The Al atoms are on the sphere of 
4.03 A radius, while the N atoms lie on the sphere of 4.18 
A. These differences in radii of Al and N skeletons are of 
same magnitude to those found for the B and N radii in 
the B24N24 cage. In contrast, the S4 cage, also shown in 
Fig. ^ satisfies the isolated six square hypothesis^ and 
does not contain any octagon. The isolated square hy- 
pothesis is similar to the better-known isolated pentagon 
rule for carbon fullerenes. The exact A1N analogue of 
carbon fullerenes do not exist. Fullerenes do not per- 
mit full alternation of aluminum and nitrogen atoms due 
to presence of pentagonal rings. Even numbered rings 
are necessary in alternate cages. Using Euler's theorem, 
it can be shown, that alternate A1N cages can be made 
closed using exactly six squares (four atom rings). The 
S4 structures is such a cage containing twenty hexagons 
and six squares. It is proposed as a fullerene-like ana- 
logue of an alternate boron nitride fullerene. 

The cage contains two octagons, eight squares, and 
twelve hexagons. It is indistinguishable from the 0cage 
when viewed along four fold axis passing through the 
center of octagon. It is basically a very short closed 
capped (4, 4) armchair aluminum nitride nanotube. As 



it is practiced in literature we use the term nanotube 
for cylindrical cage structures. Incidentally, the caps of 
(4, 4) armchair A1N NT are the hemispherical halves of 
round octahedral 0cage. 39 For boron nitride, the 5s and 
S4 cages are energetically favored over the 0cage, al- 
though the energy differences were quite small In 
the present AI24N24 case, similar trend is observed and 
the energy differences are further diminished. The S4 and 
S$ cages are energetically nearly degenerate with 0cage 
being higher by 0.1 eV. 

TABLE II: The electronic structure, and symmetry of HOMO 
and LUMO of the A1N cages and nanotubes. 



System Electronic structure HOMO LUMO 



A124N24 


O 


10ai 10a 2 30ti 30t 2 20e 


e 


Ol 


A124N24 


s 4 


60a 606 120e 


e 


a 


A124N24 


Ss 


30a 306 60ei 60e 2 60e 3 


a 


a 


A1 2 8N 28 




38a 9 32a„ 386 9 326 u 64e 9 76e u 


a u 


a g 


A128N28 


T 


26a 44e 70t 


t 


a 


A132N32 


s 8 


40a 40b 80ei 80e 2 80e 3 


b 


a 


A1 3 6N 3 6 


T d 


24ai 6a 2 30e 36ti 54t 2 


t 2 


ai 


A1 48 N48 


Ss 


60a 606 120ei 120e 2 120e 3 


a 


a 


AlgeNge-I 


O 


40ai 40a 2 80e 120ii 120t 2 


a 2 


ai 


A1 96 N96-II 


O 


40ai 40a 2 80e 120ii 120t 2 


(12 


ai 


The S s 


Al 24 


N24 cage can be extended alo 


ng the four 



fold axis by inserting a ring of alternate A1N pairs to ob- 
tain AI28N28 nanotubes with C^h symmetry. The C±h 
A1 2 8N 2 8 tube can also be generated from the 0AI24N24 
cage by cutting former into two halves, orienting them 
along the C4 axis and then inserting a ring of eight al- 
ternate A1N atoms, followed by a rotation of one half by 
forty-five degrees. Likewise, the C4/1 AI28N28 tube can be 
extended by adding one and five rings of alternate A1N 
atoms to obtain AI32N32 and AI48N48 nanotubes, respec- 
tively. Both nanotubes have Ss symmetry. The addition 
of successive alternate rings will lead to the more famil- 
iar infinite (4,4) A1N nanotube. The optimized structures 
of the capped nanotubes thus derived are shown in Fig. 
02 The tubes have diameter of 6.68 A. The largest tube 
A1 48 N48 has length of 16.5 A. The A1N bond lengths in 
tubes in general vary from 1.74 A to 1.80 A. The A1N 
bond that shares hexagonal and octagonal rings is the 



FIG. 3: (Color online) Finite A1N nanotubes: AI24N24 , 
A1 28 N28 , AI32N32 , and A1 48 N 4 8 



shortest (1.74 A) while the largest one (1.80 A) is shared 
by octagonal and square rings. The inner hexagonal rings 
have A1N bond distances in the rang e 1.76-1.79 A. 

The AI28N28 cage of T symmetry and the tetrahedral 
Al36N 36 cage are another two cages that satisfy the iso- 
lated six square hypothesis, and thus prior to the experi- 
mental boron nitride workii were considered most favor- 
able. The former contains twenty four hexagons while 
the later has thirty-two hexagons. The boron nitride 
counterpart of both these structures have been proposed 
as candidate structures for abundant boron nitride clus- 
ters detected in mass spectrum. The A1 2 8N 2 8 cage of T 
symmetry is energetically more favorable than the Cih 
AI28N28 nanotube. The former has six four-fold rings 



FIG. 4: (Color online) Octahedral AI96N96-I cage as seen 
from the four- fold (top), three- fold (middle), and two- fold 
(bottom) axes, respectively. 

that isolated by hexagonal rings while the latter has more 
defects eight squares and two octagons. The extra de- 
fects in the C\h make it energetically less favorable. Our 
calculations show that the A1N bond distances in the Td 
cage are largest (1.80 A) for the bonds that belong to the 
square ring. The bonds between a square and hexagon 
are 1.76 A in length, and those between hexagonal rings 
have length of 1.79 A. 

The Al9 6 N 96 cluster is a larger cage that preserves the 
symmetry of the AI24N24 by surrounding all square and 
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FIG. 5: (Color online) Octahedral AlgeNge-II cage as seen 
from the four- fold (top), three- fold (middle), and two- fold 
(bottom) axes, respectively. 

octagonal defects by rings of hexagons like leapfrogging 
in the fullereneSj^i which also leads to a cluster four times 
as big. Leapfrog fullerenes necessarily have closed elec- 
tronic shells. All III-V fullerene-like clusters have closed- 
shells due to large HOMO-LUMO gaps. The A1 96 N 96 
cage contains eighteen defects: twelve squares and six oc- 
tagons. The optimization of AI96N96 , starting from the 
ionic configuration of the B 96 Ng 6 cage scaled by roughly 
the ratio of A1N and BN average bond distances, results 



in the fullerene-like AlggNgg-I cage as shown in Fig. 0] 
It is evident from the figure that squares stick out. The 
cage has an average radius of 7.5 A and squares pro- 
trude out by roughly 1 A from the average radius. On 
the other hand, if the optimization of Alg 6 Ng 6 is started 
from the unsealed BggNg6 cage, then one obtains a very 
different 0structure. The resultant structure, called two- 
shell AlggNg6-II cage hereafter, is shown in Fig. [3] It is 
very different from the AlggNggT cage in Fig. 01 and is 
basically surface growth on the octahedral AI24N24 cage. 
The AI24N24 core of this structure is evident in Fig. 03 
The atoms on the surface of inner core have four-fold 
coordination. The two-shell cage is lower in energy by 
0.66 Hartrec than the fullerene-like cage despite having 
many more defects than the fullerene-like AlggNgg-I cage. 
From this result it appears that larger spherical clusters 
or quantum dots of A1N would prefer a solid-like cages 
rather than hollow fullerene like cages. On the other 
hand, earlier study have found larger AIqqNqq cage to be 
energetically and vibrationally stabled It would therefore 
be interesting to investigate relative stability of larger 
clusters containing fullerene-like and solid-like spherical 
quantum dots. We also examined the two-shell cage as 
a candidate structure of BggNg6 . Our calculation indi- 
cate that the cluster is not stable and fragments upon 
optimization. 

The electronic structures of A1N cages are given in 
Table |n] while the binding energies per A1N pair, the 
HOMO-LUMO gaps, the ionization potential and the 
electron affinity of these cages are given in Table IIIII 
The ionization potential (electron affinity) is calculated 
from the difference in the self-consistent solutions of neu- 
tral A1N cage and its cation(anion). All structures are 
energetically stable with binding energies of order of 10- 
11 eV per A1N pair. It is well known that the HOMO- 
LUMO gaps obtained from the density functional models 
including the present one underestimate true band gap. 42 
The so called A SCF provides good approximation for 
the ionization potential and the electron affinity of the 
system4& We determine the band gap as a difference of 
ionization potential and electron affinity computed by the 
A SCF method. The time dependent density functional 
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TABLE III: The calculated values of binding energy per 
A1N pair (BE), the energy gap between the highest occu- 
pied molecular orbital and the lowest unoccupied molecular 
orbital, the vertical ionization potential (VIP), the electron 
affinity (VEA), and the energy gap obtained from the ASCF 
calculation for the optimized A1N cages. These calculations 
are spin-polarized. All energies are in eV. 



Symmetry BE GAP VIP VEA A SCF 



A124N24 




O 


10.24 


2.97 


7.05 


1.46 


5.59 


A124N24 




Si 


10.34 


2.47 


6.84 


1.72 


5.12 


A124N24 




S 8 


10.34 


2.63 


6.79 


1.58 


5.21 


A128N2S 




Cah 


10.42 


2.74 


6.81 


1.59 


5.22 


A1 28 N28 




T 


10.45 


2.67 


6.84 


1.69 


5.15 


A132N32 




Sa 


10.49 


2.79 


6.77 


1.61 


5.16 


A1 36 N 3 6 




T d 


10.54 


2.70 


6.73 


1.76 


4.95 


AI48N48 




s d 


11.09 


2.81 


6.56 


1.76 


4.8 


A1 96 N 9 6- 


-I 


O 


10.56 


2.63 


6.22 


2.04 


4.12 


AlgeNge- 


-II 


O 


10.75 


1.81 


6.03 


2.53 


3.50 



theory (TDDFT) or the GW approximation provide cor- 
rections to the HOMO-LUMO gaps and are more suitable 
for calculations of the band gaps. We have not yet imple- 
mented these techniques. However, the present method- 
ology of using fitting basis sets can also make TDDFT 
calculations efficient^ The A SCF calculated gaps are 
given in Table II I II In case of AI24N24 although the O 
cage is energetically less favorable its HOMO-LUMO gap 
and A SCF gap are larger than those of the S4 and Ss 
cages. Amongst the A1N cages that have only six isolated 
squares, the binding energy systematically increases for 
S 8 AI24N24, T A1 28 N28, and A1 36 N 36 cages. This trend is 
similar to that observed in case of boron nitride cages^S 
and carbon fullerenes4^ The squares stick out too far 
in the M96N96 (M= B, Al) cage for that cage to be fa- 
vored and its hemisphere to be an extremely favorable 
cap of the (8,8) nanotubes. Both the AlggNgg cages have 
smaller binding energy than the AI48N48 nanotube, per- 



haps due to larger number of defects (octagonal and four- 
fold rings). The O AI24N24 and AL30N60 cages were also 
studied by Change etali^i They find that former has bind- 
ing energy of 4.72eV/atom while the latter has binding 
energy of 4.76 eV/atom. The increase in the binding en- 
ergy from octahedral A1 2 4N 2 4 cage to AL30N60 cage is 0.05 
eV/atom. Our calculations show binding energy gain of 
0.4 eV/atom with size increase from octahedral AI24N24 
to octahedral AlggNgg. It is therefore likely that the oc- 
tahedral cages studied in this work are more stable than 
the icosahcdral AlgoNgo cage. In comparison with their 
boron nitride counterparts^S^^I the A1N cages are ener- 
getically less stable, have lower ionization potentials and 
smaller band gaps. This feature is similar to the bulk 
phase of these materials. The band gap of solid A1N is 
smaller than that of BN solids. 

To summarize, fullerene-like cage structures and finite 
(4, 4) nanotubes are studied by the density functional cal- 
culations using polarized Gaussian basis sets of double 
zeta quality. The binding energy, electron affinity, ion- 
ization potential, the HOMO-LUMO gap and the ASCF 
gap are calculated for the optimized A1N cages. Calcu- 
lations show that all A1N cages are energetically stable, 
with band gap of order of 5 eV. For the larger AlggNgg 
cluster, the two-shell cage with an interior AI24N24 0cage 
is energetically favorable over the fullerene-like cage. The 
binding energies and band gaps are smaller than their BN 
counterparts. We hope that the present study will aid the 
experimental search for the A1N cages. 

The Office of Naval Research, directly and through 
the Naval Research Laboratory, and the Department of 
Defense's High Performance Computing Modernization 
Program, through the Common High Performance Com- 
puting Software Support Initiative Project MBD-5, sup- 
ported this work. 
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